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Abstract 

The Euclidean path integral method is applied to a quantum tunneling model which accounts 
for finite size (L) effects. The general solution of the Euler Lagrange equation for the double well 
potential is found in terms of Jacobi elliptic functions. The antiperiodic instanton interpolates 
between the potential minima at any finite L inside the quantum regime and generalizes the well 
known (anti)kink solution of the infinite size case. The derivation of the functional determinant, 
stemming from the quantum fluctuation contribution, is given in detail. The explicit formula for 
the finite size semiclassical path integral is presented. 

PACS numbers: 03.65.Sq Semiclassical theories and applications, 03.65.-w Quantum mechanics, 31.15.xk 
Path Integral Methods 
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1. Introduction 



As a fundamental manifestation of quantum mechanics the tunneling through a poten- 
tial barrier plays a central role in many areas of the physical sciences. Being intrinsically 
nonlinear phenomena, quantum tunneling problems cannot be attacked by standard pertur- 
bative techniques while semiclassical methods are known to provide an adequate conceptual 
framework [l, 2]. 

To be specific, I consider a particle of mass M moving in a bistable symmetric double 
well potential with oscillation frequency u: 

rr/ \ , 5 4 5 4 ... 

V(x) = —x 2 + -x 4 + -a 4 (1) 

The minima are located at x — ±a and the positive quartic force constant 5 (in units 
eVA~ 4 ) is related to the other potential parameters by 5 = Muo 2 /a 2 . The latter term in 
Eq. (CO) ensures that the potential is positive defined and V(x = ±a) = 0. As it will be 
shown in detail below, the minima are asymptotically connected by the classical (anti)kink 
solution, a charge conserving domain wall whose energy is inversely proportional to 5 for a 
given vibrational frequency inside the well. Thus, even small quartic nonlinearities induce 
large classical energies which, in turn, enhance the effects of the quantum fluctuations in 
the overall probability amplitude for a particle to move from one minimum to the next. 
This explains why, on general grounds, perturbative methods fail to describe the tunneling 
physics. 

In general, 4 models also host potentials having a quartic 5 < [3] and positive sign of 
the quadratic term in Eq. (JT]). In this case the mid-point of the potential valley would be a 
classically and locally stable configuration, a metastable state from which the particle would 
tend to escape via quantum tunneling. In such state, known in many branches of physics 
as "false vacuum" [4J, the system has a finite lifetime due to the presence of a negative 
eigenvalue in the quantum fluctuation spectrum which governs the decay into the "true 
vacuum". The classical solution of the Euler-Lagrange equation for metastable potentials 
is a time-reversal invariant bounce carrying zero charge. Also this object, whose physical 
properties generally differ from the (anti)kink, has an energy inversely proportional to the 
anharmonic parameter 5 which makes perturbative analysis of metastability unfeasible. 
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It has been noted since long that the 5=0 singularity in the thermodynamic limit of 
the partition function for the anharmonic potential is formally similar to the condensation 
point in the droplet model for phase nucleations |3|. The nature of such singularity has 
been widely studied in the past also in connection with fundamental investigations of the 
anharmonic oscillator and its energy levels 0,0,0. 

Tunneling problems in quantum field theory have been treated by semiclassical methods 
9| which extend the WKB approximation [10| and determine the spectrum of quadratic 
quantum fluctuations around the classical background. While such methods had been first 



applied to infinite size systems, generalizations to finite sizes |U[ have been later on devel- 
oped to obtain correlation functions and spectra of 4 potentials The Casimir effect in 
4> A theories is another remarkable example in which the properties of the system depend on 
the size of the cavity 
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Finite size effects have become a focus of research also on classical mechanical systems 
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16] in which spatio-temporal noise of classical origin induces escape from a metastable 



state. Examples of current technological interest are the switching rate in the magnetization 
of micromagnets and the finite lifetime of metallic nanowires shaped by cylinders whose 
radius undergoes fluctuations of thermal origin 18 1. Such systems exhibit a structure of 
thermal activation regimes which crucially depend on the system size (L) and the boundary 
conditions. The transition in the classical activation rate versus L is formally identical to 
the crossover from the classical activation to the quantum tunneling regime once L is made 
proportional to the inverse temperature of the system. 

In general, the finite size formalism offers the bridge to derive the thermodynamics of 
the system in the spirit of the thermodynamic Bethe Ansatz [19] : this is accomplished by a 
Wick rotation which transforms the real time t into the imaginary time r of the Matsubara 
formalism, t — > —it, where r varies in a periodic spatial box of finite size L = h(3 with 
(3 = (KbT)^ 1 , Kb being the Boltzman constant and T the temperature. Accordingly the 
infinite size system maps onto the zero temperature limit. 



With these premises, the path integral method in the Euclidean version |20l . l2l| emerges 
as a natural tool to deal with finite size/temperature tunneling phenomena in semiclassi- 
cal approximation. The quantum probability amplitude is given in terms of a functional 
integral over all histories of the system depicted by those r— dependent paths which ful- 
fill (anti)periodic boundary conditions over a size L. In the semiclassical treatment, the 
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functional integral reduces to Gaussian path integrals over quadratic quantum fluctuations 
around the classical solution of the Euler-Lagrange equation. The fluctuations contribution 
to the amplitude can be evaluate d by a number of techniques in the context of the func- 



tional determinants theory 



22 



23 



24] that implements the Euclidean path integral method. 



While the latter is well known for the infinite size system (both for bistable and metastable 

n 

4 models) [25j], much less known is the mathematics and the associated physics for the 
finite size case. A contribution to the research field comes from this paper which analyses 
the finite size Euclidean path integral for the bistable 4 model proposing a new family of 
classical backgrounds. I emphasize that the model is intrinsecally non dissipative as there 
is no coupling to the heat bath while the temperature, as stated above, provides a measure 
of the system size. The work is organized as follows. Section 2 contains the solution of 
the Euler-Lagrange equation, the finite time instanton which minimizes the Euclidean ac- 
tion. This family of paths, expressed in terms of Jacobi elliptic functions, generalizes the 
(anti)kink peculiar of the infinite time theory and interpolates between the minima of the 
potential in Eq. (CQ). Section 3 presents the generalities of the semiclassical path integral 
formalism posing the problem of the evaluation of the quantum fluctuations functional de- 
terminant. The latter is derived in terms of elliptic integrals and analysed in Section 4. 
Some final remarks are made in Section 5. 



2. Classical Background 

Consider the potential of Eq. (TJ|). The real time formalism does not allow classical 
motion in a double well potential but, switching to the imaginary time formalism t — > —ir, 
the potential is reversed upside down and the classical orbit x c i(t) is obtained by solving 
the classical equation of motion: 

Mx d (r) = V'(x d (r)) (2) 
where V means derivative with respect to x c \. Integrating Eq. (J2j) , one gets 

^W-^,(r))= E (3) 
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with integration constant E. From Eq. ([3]) one easily derives the equation of motion in 
integral form: 



_ ± [M_ M t) dx 
T ° V 2 J Xcl{To) ^ E + V(x) [ J 

The zero energy configuration is consistent with a physical picture in which the particle 
starts on the top of a hill at x c i(t = — oo) = —a with zero velocity, crosses the valley and 
reaches the top of the adjacent hill at x c i{t = +00) = +a with zero velocity. r represents 
the time at which the bottom of the valley (of the reversed potential) is crossed and, by 
virtue of the r— translational invariance, it can be placed everywhere along the imaginary 
time axis. This freedom has a price to be paid, that is a zero mode whose eigenvalue causes 
a divergence in the path integral as discussed below. 

Imposing x c i(tq) = 0, one gets from Eq. (j4j) with E — 0, the well known (anti)kink 



x c i(t) = ±atanh[co>(r — t )/\/2] (5) 



Since the transition happens in a short time, almost instantaneous 



y, the (anti)kink so- 



lutions are also named (anti)instantons in quantum field theory [2f], [27]. There are however 
other solutions of Eq. (jl|) associated to finite values of E and corresponding to the physical 
requirement that the classical path has to connect the potential minima at finite times. Thus, 
the finite size problem needs to be solved with antiperiodic boundary conditions (APBC) 
such that x c i(r = ±L/2) = ±a. Define 



Xd{r) - 



V2 a 
2E 
2 ' 5a 1 



1 2£ 

+ 7-I ( 6 ) 



Then Eq. (jl]) transforms into 



1 d X m 

r " r o = ±- / t A „ . 7 = (7) 



W Xci( ro) vV - X 2 + k/2 
Solutions of Eq. ([7]) can be found in three ranges of k: 



1) k < 0. Xdij) is unbounded and always larger than a, hence it does not fulfill our physical 
requirements. 

2) < k < 1/2. In this range and, for \x\ < \J ^^Y^ 1 ' there is a bounded, continuous and 
antiperiodic family of solutions given by 



x d (r) = ±clF(k) ■ sn [ V 71 + J^^ ^^( r _ 



2 « x v2k 



r = ; — ; F(« 



(8) 

This solution has often been proposed in the finite size studies on the double well potential 



28l . |29| . l30f | . Around this background one can straightforwardly determine the low lying 
excitations in the fluctuation spectrum by solving a standard Lame equation [31]. Note 
however that F(k) in Eq. (jSJ) is always smaller than one. As the Jacobi elliptic sn-function 
is < 1, the classical paths in Eq. (jSJ) undergo periodic oscillations between the slopes of the 
hills in the reversed potential but never join the potential extrema (unless the trivial case 
k — 1/2 is recovered). This is physically due to the fact that k < 1/2 implies E < (see 
Eq. ©) but negative energies are not consistent with finite path velocities on the r— range 
boundaries. In fact, from Eq. ([3]) it is easily seen that the relation 



x d {r = ±L/2) = ^2E/M (9) 

holds on the r— boundaries. This observation alone proves that, unlike it has been often 
assumed in the past, paths with E < cannot describe the physics of the finite size double 
well potential. Therefore Eq. (jSJ) is useless to our aim. 

3) 1/2 < k. This turns out to be the physically interesting case since in this range the 
classical motion has E > 0. The solution of Eq. ([7]) is presented in detail. 
The integral in the r.h.s. of Eq. (j7j) yields: 



6 



dx 



Vx 4 ~X 2 + k/2 
2 air-tan ( W — x 



F(9,s) 



5<" 



(10) 



where -F(#, s) is the elliptic integral of the first kind with amplitude 9 and modulus s: 



Jo y 1 — s z sin a 

whose quarter-period is given by the first complete elliptic integral K(s) = F(n/2,s) 
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321 ] . In terms of F(6,s), one defines the Jacobi elliptic functions sn-amplitude, cn- 



amplitude, dn-amplitude : 



sn[F(6,s)) = sin9 
cn[F{6,s)} = cos# 

dn[F(6,s)] = Vl - s 2 sm 2 6 (12) 



From Eqs. (JTj), (TTDT) one gets 



and using the relations: 



±\/ -2uj(t - r ) = F(B,s) (13) 



2tan(0/2) 
tan# - v 1 ' 



1 - tan 2 (#/2) 
m[F(9,8)] 
cn[F(6,s)} 



tane = f!iM (14) 



Eq. (TTB"j) transforms into: 



9 4 /~.-v 

sn(2o7, s) z ' V 
cn(2ru, s) -, _ /? 2 



07= ± 4 A /^(r-r ) (15) 
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The condition 4 y fxd( r ) 7^ =tl imposed by Eq. ([To']) is equivalent to x c i(r)/a ^ ±(2s 2 — 
l) -1 / 2 which is always fulfilled in the physical range of our interest. The l.h.s. in the first of 
Eq. (I15p can be rewritten using the double arguments relations: 



2sn(zu, s)cn(zu, s)dn(zu, s 
sn(2zu, s) = — 



1 — s 2 sn 4 (zu, s) 
cn 2 (w, s) — sn 2 (zu, s)dn 2 (zu, s) 



1 — s 2 sn 4 (zu, s) 

(16) 



cn(2zu, s) 



and, through the definition in Eq. ([6]), from Eqs. ([15]) . ([16]) I derive the two solutions for 
the classical equation of motion in the finite size model: 



(l) / \ 4 nr sn(zu,s)dn(zu,s) 



cnizu, s) 



(2), x 4 nr cn(m,s 
ci v ) = - v2«a 



sn(w, s)dn{w, s) 

(17) 

The antiperiodic paths of Eq. (fTTj) are plotted in Fig. [1] and Fig. [2] respectively, for 
a choice of potential parameters ui, M, a and energy E which fulfill the criteria for the 
validity of the semiclassical approximation. Aside from the prefactor 4 v^2/ta, the two families 
of classical paths are antireciprocal of each other. Thus the zeros of x^\t), occuring at 
zu = ±2nK(s) with integer n would be divergencies for x^ (t). Vice- versa, the zeros of 
x d ( r ) a ^ w = =t(2^ + ^-)K(s) are singular points for x^\t). 

In fact such singularities, that remind us of the caustic points in the classical orbit of the 
harmonic oscillator [33], happen to be outside the range physically relevant to our purposes. 
I clarify this point. 

First, note that the solutions in Eq. (jTTj) have a period 4K(s) thus, in principle, they 
should be discussed in the range w G [— 2K(s), 2K(s)) since the amplitude 6 of the elliptic 
integral varies as 6 G [— 7r, 7r]. Second, focus on x^\t) and apply the APBC, x^\r = 
±L/2) = ±a, emboding the physics of the finite size problem. Setting tq = without loss 
of generality, the first of Eq. ( fTTj) yields: 
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FIG. 1: (Color online) Classical path x\, (r) in Eq. (|17p . for a choice of energy i£ (in meV) and 
potential parameters: a (in A), fiiu (in meV) and particle mass M in units of the electron mass. 
The center of the (anti)inst 
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FIG. 2: (Color online) Classical path x , (t) in Eq. (|17p . for the same parameters as Fig. [TJ 



4 / — sn(cc7 , s)dn(w , s) 

1= V2k t r 

cn(zu*, s) 



Representing the elliptic functions in Fourier series [311 ] . I calculate Eq. ffTTj) and obtain 
the values ro* fulfilling the APBC. Precisely, the computational method sets a value of E and 
determines the corresponding w* such as Eq. f fT8|) holds. In this way the program calculates 
the period L = L(E) and establishes the one to one correspondence between E and u* at 
which the path interpolates between the minima. The ratio w*/K(s) turns out to be always 
smaller than one, accordingly the singularities at w — ±K(s) are physically avoided. In 
other words, the path (r) always connects the potential minima well before encountering 
the singular points. Thus, taking the center of motion at x c /(r = 0) = 0, x^\t) properly 



describes a family of finite size (anti)instantons which are bounded, continuous and odd in 
the range r G [-L/2,L/2). 

An analogous procedure with APBC applied to x^\t) in Fig. pleads to determine the 
boundaries ±m*/K(s). As r may be located everywhere along the r— axis one may define 
in principle an even range m/K(s) with no singularities also for (t). However, setting 
tq = 0, the choice of x^\t) happens to be more convenient since the latter is well behaved 
for r G [-£/2,L/2]. 

The remarkable fact is that the (u*,E) correspondences obtained by the two families 
of classical backgrounds are essentially identical. This means that the two independent 
solutions of the finite size classical equations of motion predict the same physics. 

Correctly, in the limit E — > (k — > 1/2), one gets from Eq. (fT7|) : 



x^\t) -> ±atanh[w(r-r )/V2] 
x d ( T ) ~~ > =Facoth[u;(r - r )/\/2] 



(19) 



The first in Eq. (fT9l) is the (anti)kink whose center of motion is set at x^\t ) = whereas 
the second is the solution of the infinite size classical equation Eq. (j2J) fulfilling the condition 

%cl ( r o) = OO- 

Further, using the relations: 



d sn(w, s) 

dm 
dcn(zu, s) 

dm 
d dn(m, s) 

dm 



cn(m, s)dn(m, s) 
—sn(m, s)dn(m, s) 
—s 2 sn(m, s)cn(m, s) 



one gets, from the first in Eq. (|17|) . the particle path velocity as: 



(20) 



Xcl{T) 



duly K 



22 sn 2 (m,s)dn 2 (m,s) 
1 — 2s sn (m, s) + - 



cn 2 (m, s) 



(21) 



Using the second of Eq. (1151) . one can easily check that Eq. (12111 does not vanish on the 
r-range boundaries consistently with the general constraint expressed by Eq. (Q. The latter, 
together with Eq. ([21]) . permits to define the one to one correspondence between E and u*. 
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3. Semiclassical Path Integral 



In the semiclassical model, the particle path x(t) splits in the classical background x c i(r) 
plus the quantum fluctuations i](t). Thus, the path integral for space-time particle propa- 
gation between the positions Xi and Xf in an imaginary time L is: 



< Xf\Xi >l= exp 



A[x cl 

h 



Dt)(t) exp 



h 



The classical action A[x c ;] is formally given by 



A[x d ] = V2M j dx dy /E + V{x d ) - E ■ L(E) 
Inserting the potential in Eq. (flj) , Eq. (l23|) transforms into: 

A[ Xcl ] = A Q J' dj,^ + (l+y)2(l-y)2 - E ■ L(E) 

V25a A 



(22) 



(23) 



(24) 



suitable for computation while, in the limit E — > 0, one obtains the well known result: 



A[x d ] 



2V2M 2 u 3 



(25) 



3 S 

Eq. (1251) represents the (anti)kink energy (in units of H) whose inverse dependence on 5 
has been pointed out in the Introduction. 

Being quadratic in the quantum fluctuations, Af\n] defines a stability equation whose 
eigenvalues represent the fluctuation spectrum 25j. Carrying out Gaussian path integrals 
over the fluctuation paths, one gets formally 



Drj(r) exp 



JAM 



tt-Det 



O 



-1/2 



d = -d 2 T +V"{x d {r))/M 



(26) 
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where K depends on the measure of integration. The heart of the matter is the evaluation 
of the functional determinant in Eq. (126]) . Three observations permit to get the clue to the 
problem. 

A) Differentiating Eq. (J2J) with respect to r, one notes that x c i(t) solves the homogeneous 
equation associated to the second order Schrodinger-like differential operator O 



6t){t) = (27) 

As a general consequence of the r-translational invariance, a Goldstone mode appears in 
the system and this is true both for the infinite and the finite size cases. Thus, the particle 
path velocity is proportional to the fluctuation which makes Af[rj] vanishing, x c i{t) oc r/ (r) 
with eigenvalue Eq = 0. 

B) As x c i(t) is monotonic within a period L, x c i(r) has no nodes hence, T]o(t) represents 
the ground state and Eq is the lowest eigenvalue in the fluctuation spectrum in agreement 
with the fact that the potential is stable. However, the zero mode eigenvalue breaks the 
Gaussian approximation and makes the path integral divergent. The trouble is overcome 
by regularizing the fluctuations determinant, extracting e 
contribution separately resorting to collective coordinates 
for the infinite and finite size theories respectively: 



rom Eq. (1261) and evaluating its 



5 



251 ] . The replacement yields, 



(eq) — > \ f L infinite size 
V zirn 



(28) 



Note incidentally that also the infinite size instantonic theory assumes, strictly speaking, 

—1/2 

that the size L along the r— axis is finite (although large) so that e oc L and the path 
integral is finite. Otherwise the zero mode singularity could not be removed. This obser- 
vation points to a formal contradiction of the infinite size theory that, in fact, stimulated 
me to develop a consistent description in which both the classical background and the zero 
mode of the quantum fluctuations can be precisely determined at any finite size. 
C) As stated below Eq. (j2T|) . x c /(r) takes the same values on the r— boundaries. More 
generally, for any two points vj\,vji such that w 2 = W\ ± 2K(s), from Eq. (1211) one gets: 
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Xdi^2) = Xcifai). As the path velocity is (proportional to) the ground state fluctuation, the 
important consequence is that the fluctuation eigenmodes obey periodic boundary conditions 
(PBC). 

Collecting the informations from AJjto C) and using a fundamental result of the theory 



of functional determinants with PBC 

as 



231 ]. I rewrite the fluctuation determinant in Eq. (l2T)j) 



Det[6] = e -Det R [6] 



Det H [0} 



(29) 



Where Det R [0] indicates that the determinant is regularized 24j after extracting the 
zero eigenvalue, rjo and rji are independent solutions of Eq. (Ej 



Vl =^m=s> (30) 

W(r]o,f]i) is their Wronskian and < f?o|??o > is the squared norm. Using the first in 
Eq. (IT7|) and Eq. fl30l) . Det R [0] can be worked out analytically. 



4. Regularized Fluctuation Determinant 

I show now in detail how to obtain the regularized Det R [0] of Eq. (|29l) on the base of 
the Wronskian construction. The two points tu 2 , W\ map onto L/2 and — L/2 along the r— 
axis according to the definition in Eq. (I15p . 

I) The r— derivatives of the two independent solutions r] and rji of Eq. (130]) read respec- 
tively: 



_d_ 

dm 



dx c i(r) 
dzu 



dw 
Ifr. 



(31) 



Then, the Wronskian W(t] ,rii) is: 
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(32) 



The Wronskian is constant and therefore can be calculated in any convenient r G 
[— L/2,L/2]. The center of the classical paths, located at r , is the best choice as 
Xd(r ) = x d (r ) = and x d (r ) = aw^fn. 

Observing that: 



dw 



U! 



Or 2^s 2 - 1/2 



_d_ 

dm 



dw 



1/2) 



-3/2 



(33) 



from Eqs. (ETJ) and ([32]), I get 



2^2 



w a 



16 



[s 2 - 1/2)- 



(34) 



II) I evaluate 7/1 (r = ±L/2) writing 7/x as 



ro = -F(e,s) 



(35) 



Using the relations: 



dw 1 dF(6,s) 
dm 4s ds 



dF(6, s) 
ds 

E(9,s)-- 



-2 1 2 
S = 1 — S 



S 2 



E(6,s) -s 2 F(9,s) 



day 1 — s 2 sin 2 a 



sin # cos 8 
'\ - s 2 sin 2 6 



(36) 
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I derive the following boundary properties: 



9(L/2) = 2arctan^j^=^j = -9(-L/2) 

F(e(L/2),s) = -F(9(-L/2),s) 
E(9(L/2),s) = -E(6(-L/2),s) 



(37) 



After setting r = 0, it follows that: 



m (L/2) = -Vi(-L/2) 
Vi(-L/2) = x c (-L/2) 



( dw 
Vdr 



dzui 
dm 



-L/2 



b c (±L/2) = 2auo^Ks 2 {\ - sn 2 {w*)) 



x 



sn [w 



1 - cos[0(L/2)] 



1 - s 2 sin 2 [0(L/2)] 



uoL 



zu 



(38) 



Ay/s 2 - 1/2 

The last in Eq. (|38|) coincides with w* in Eq. (ITS1) . Note that 771 is a solution of Eq. (J2j 
but, unlike rj , it does not represent a quantum fluctuation component. This explains why r)\ 
is odd on the boundaries whereas rj fulfills PBC. Thus, Det R [0] in Eq. (1211 can be rewritten 
as: 



W^l^O; Vl) ^ V T ' Om t=L/2 

(39) 

The last factor in Eq. (139]) is computed via Eqs. (|36l) . (137|) and turns out to be positive 
thus ensuring the correct (negative) sign to Det R [0] (due to W(r] ,r]i) < 0). It is well 
known from the general theory that only ratios of functional determinants have physical 
meaning both in value and sign. In fact being a product of an infinite number of eigenvalues 
with modulus larger than one, each determinant is separately divergent in the L — > 00 limit. 
To get a finite ratio, I normalize Det R [0] over the harmonic oscillator determinant Det[h] 
which, for periodic boundary conditions, reads: 
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Det[h] = -4sinh 2 (c;L/2) 
h = -d 2 +u 2 



(40) 



As Det[h] and Det R [0] have the same exponential divergence, their ratio is finite in the 
L — > oo limit. Moreover it is positive, consistently with the fact that the contribution by 
the zero mode eigenvalue (extracted from Det[0}) is positive. Therefore also Det[0]/Det[h] 
has to be positive since the system is stable and there are no lower eigenvalues below the 
zero mode eigenvalue. 

Ill) f]o is proportional to the ground state ortonormal component in the series expansion 
for the path fluctuation. I define fj = Nr] such that, in Eq. fl39l) . < r] \r]Q >= N~ 2 . Then, 
using Eq. (I2TI) . I derive the squared norm in the form 



AT 



G I dw 
'o 

oja 2 



22 sn 2 (w,s)dn 2 (zu,s) 

1 -2s sn (w, s) H —. r 

en (w, s) 



G 2(s 2 - 1/2) 3 / 2 (41) 
suitable for computation. Alternatively, one may proceed through the definition of the 
classical action in Eq. ff23l) and get: 



N- 2 = -(A[x d ] + E.L{E)) (42) 

In any case < r^ol^o > has the dimension of [uia 2 ] hence, from Eqs. (1531) . (154]) . (j4Tp . it 
follows that Det R [0] is proportional to lo~ 2 and correctly it does not carry any dependence 
on a. This completes the analysis since all ingredients required to compute Eq. (I39p are now 
known. 

Taking the ratio between the determinants in Eq. fl39|) and Eq. fj40l . in the L — > oo limit, 
I get 



JL (43) 
Det[h\ 12u 2 V 

thus recovering the well known value of the infinite size instantonic approach. This proves 

the correctness of the analytical procedure. 
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Then, the semiclassical path integral in Eq. (1221) for one (anti)instanton takes the final 
expression: 



exp 



< Xf\Xi > L 

Q{L) = 
n QF (L) 

which can be evaluated using Eq. 



M 



27rhL 2sinh(a;L/2) 
A[xd] 



tt{L)L 



h 



MN-- 
2nH 



n QF (L) 

Det R [6] 

, (HOD, (S2). Note that: i) The harmonic 

M 
2-kKL 



determinant is, in turn, normalized over the free particle determinant Det[— d%] = 
which incorporates the constant K. ii) Vfi F (L) accounts for the quantum fluctuation effects 
in the finite size tunneling problem, in) Q(L) is the overall tunneling frequency which 
removes the twofold degeneracy of the double well potential in a finite domain. From the 
path integral in Eq. (|44p one can extract the physical properties of the closed quantum 
system. 



5. Final Remarks 

This work presents the mathematical description of the finite size semiclassical theory 
for the bistable <ft 4 potential based on the Euclidean path integral formalism. In particular 
the two key ingredients of the model, namely classical equation of motion and quantum 
fluctuation spectrum, are analysed in detail. I have studied the classical equation of motion 
and obtained the family of classical paths which interpolates between the potential minima 
in the finite size system: such paths are necessarily associated to positive classical energies 
and fulfill antiperiodic boundary conditions. Around such non trivial background one may 
derive the mass spectrum in 4 theories 34| by solving a stability equation more general 
than the standard Lame equation associated to Eq. ([8]) 3^]. 

Due to the time translational invariance, the spectrum of the quadratic quantum fluc- 
tuations contains a zero mode which breaks the Gaussian approximation and makes the 
path integral divergent. The result of the regularization procedure for the finite size case is 
presented. Using the theory of the functional determinants, I have derived the regularized 
fluctuations determinant giving the full set of equation required to compute it. Finally, I 
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have evaluated the contribution to the path integral originating from the quantum fluctua- 
tions and obtained an explicit expression for the tunneling energy in the finite size bistable 
(f) A model. This study permits to compute, for specific choices of potential parameters, the 
finite size induced renormalization of the tunneling energy with respect to the predictions 
of the infinite size instantonic approach. In view of the mapping to the temperature axis 
allowed by the Matsubara formalism, the obtained results may be well applied to calculate, 
for instance, the thermodynamical properties of systems described by Ginzburg-Landau field 
theories [36|, |37j, the finite T splitting energy in two level systems arising in amorphous com- 



pounds 38j, [39] and materials with local structural instabilities 



40 



41 



421 ] as well as in solids 



showing macroscopic quantum tunneling of magnetization 43]. 

Along similar semiclassical patterns one can study a metastable <ft 4 model in finite size. 
For the latter one has to build the generalized bounce solution (whose properties however 
differ very much from those of the instanton) and analyse the quantum fluctuation spectrum 
with particular care to the ground state negative eigenvalue which causes metastability. The 
softening of the low lying fluctuation modes together with the decay rate for a finite size 
system will be the subject of a next work. 
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